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The  fundamental  chemistry  and  physics  of  munitions  under  extreme  conditions 

ABSTRACT 

We  developed  the  ReaxFF  first-principles  based  reactive  molecular  dynamics  (RMD)  modeling  approach  to  determine  the  nanoscale 
phenomena  underlying  shock  detonation  processes  of  energetic  materials  (EM).  Using  the  ReaxFF  approach,  we  proposed  Compressive 
Shear  Reactive  Dynamics  (CS-RD)  simulation  methodology  to  predict  sensitivity  of  explosive  crystals  under  combined  shock  and  shear 
load.  We  also  developed  empirical  van  der  Waals  correction  to  Density  Functional  Theory  for  calculating  accurate  equation  of  states  (EOS) 
of  EM.  We 

implemented  ReaxFF  in  parallel  multiprocessor  software  to  carry  out  large-scale  simulations  of  initiation  chemistry  in  homogeneous  and 
heterogeneous  HE  under  mechanical  shock  and  shear  on  supercomputers.  We  discovered  that  sensitivity  is  dominated  by  a  combination  of 
shear  and  compression,  with  the  rate  of  decomposition  and  temperature  increase  correlating  with  the 

experimental  differences  in  sensitivity.  The  second  major  focus  was  on  the  development  of  multiscale  modeling  of  HE  detonation  using 
novel  finite  elements  method  with  explicit  generation  of  slip  lines  in  the  subgrain  microstructure  and  inclusion  of  thermochemical 
constitutive  parameters  obtained  from  RD  modeling  to  predict  the  hot  spot  formation  and  reaction  initiation  at  the  subgrain  scale  in 
polycrystalline  explosives.  The  methodologies  were  successfully  tested  and  validated  by  computational  prediction  of  anisotropic  sensitivity 
of  PETN,  as  well  as  the  formation  of  hot  spots  and  chemical  initiation  in  poly  crystalline  PETN. 
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Objective 

Develop  modeling  and  simulation  technologies  validated  by  experiment  and  sufficient  to  predict 
the  chemical  and  physical  phenomena  involved  in  the  energy  initiation  and  release  processes  of 
high  explosives  (HE)  under  extreme  conditions  of  mechanical  and  thermal  energy  deposition. 
Our  goal  is  to  develop  methodologies  for  large-scale  reactive  dynamics  (RD)  simulations  to 
determine  the  atomistic  mechanisms  relevant  to  detonation  and  sensitivity.  Elucidating  the 
interplay  between  the  reaction  chemistry  and  the  shock  dynamics  is  important  for  understanding 
the  shock-to-detonation  transition  when  reaction  products  may  affect  the  mechanical  strains  and 
stresses  in  reacting  material  and  resulting  shear  /  heat  flows,  and  vice  versa.  One  of  the  critical 
issues  is  how  the  strength  properties  and  failure  mechanisms  (shear  banding,  slipping,  fracture, 
etc)  of  the  material  influence  the  dissipation  of  the  mechanical  energy  of  impact.  This  requires  a 
detailed  microscopic  description  of  the  anisotropic  nonlinear  response  of  the  material  under  large 
deformations  to  determine  preferred  slip  planes  and  shear  flows  in  HE  micrograins  under  stresses 
generated  by  the  interference  of  shock  waves  and  the  material  microstructure 

Approach 

While  the  propagation  of  the  steady  detonation  front  is  accurately  predicted  by  continuum 
mechanics  models,  how  this  detonation  front  forms  due  to  an  initial  mechanical  or  thermal 
excitation  is  not.  This  is  because  the  initiation  process  depends  strongly  on  sub-continuum 
properties,  such  as  slip  and  shear  bands  formation,  grain  orientation,  intergranular  interfaces, 
voids,  inclusions,  and  crystal  morphology.  In  attempt  to  elaborate  the  mechanisms  of  energy 
transfer,  excitation  and  decomposition  of  energetic  molecules  we  suggest  to  use  first-principles 
based  atomistic  simulations  capable  to  capture  basic  physical  and  chemical  processes  dominating 
initiation  of  HE  and  determine  critical  properties  affecting  their  sensitivity. 

The  hot  spots  formation  in  heterogeneous  materials  usually  occurs  due  to  the  interaction  of 
shocks  with  material  heterogeneities  such  as  grain  boundaries,  defects,  voids  and  cracks.  In  case 
of  detonation  in  homogeneous  materials,  the  initiation  mechanisms  may  be  related  to  the 
heterogeneity  of  crystal  plasticity  at  the  subgrain  scale  with  high  localization  of  strain,  stress,  and 
temperature  extremes.  Conventional  mean-field  models  convey  no  information  about  these  local 
extremes  and,  therefore,  are  insufficient  for  present  purposes. 

In  this  work,  we  have  developed  a  multiscale  model  that  explicitly  accounts  for  three  scales:  (i) 
the  polycrystalline  structure  at  the  macroscale,  (ii)  single  crystal  plasticity  (including  subgrain 


microstructure  formation)  at  the  mesoscale  and  (iii)  chemical  kinetics  at  the  molecular  scale.  The 
centerpiece  of  the  model  is  an  explicit  construction  that:  i)  gives  the  effective  or  macroscopic 
behavior  of  plastically  deforming  crystals  with  microstructure,  and  ii)  enables  the  reconstruction 
of  optimal  microstructures  from  the  computed  macroscopic  averages.  An  intrinsic  feature  of  the 
optimal  deformation  microstructures  is  the  presence  of  highly  localized  regions  of  plastic 
deformation,  or  slip  lines  (see  Fig.  2.1).  Temperatures,  strain-rates  and  pressures  in  these  slip 
lines  rise  well  in  excess  of  the  average  or  macroscopic  values.  Slip  lines  thus  provide  a  plentiful 
supply  of  likely  initiation  sites,  or  hot-spots,  in  defect-free  crystals. 

We  develop  and  use: 

1.  The  reactive  force  field  ReaxFF  framework  for  describing  the  complex  chemical  and  energy 
releasing  phenomena  in  play  during  the  initial  steps  of  the  interactions  between  the  threat  and  the 
HE  for  both  shock  impact  and  thermal  impact  over  the  full  range  of  pressures  and  temperatures. 

2.  Empirical  vdW-correction  of  the  density  functional  theory  (DFT)  for  accurate  calculations  of 
constitutive  properties  and  equation  of  state  (EOS)  of  HE  over  a  wide  interval  of  strains  to 
achieve  accurate  description  of  their  response  under  shock  /  thermal  loading. 

3.  Compressive  Shear  Reactive  Dynamics  (CSRD)  to  predict  the  anisotropy  of  shock  sensitivity 
on  HE  single  crystals. 

4.  Large-scale  parallel  codes  for  RD  modeling  of  initiation  chemistry  in  homogeneous  and 
heterogeneous  HE  under  thermal  and  mechanical  shocks  and  shear. 

5.  Multiscale  model  that  explicitly  accounts  for  three  scales:  (i)  the  poly  crystalline  structure  at 
the  macroscale,  (ii)  single  crystal  plasticity  (including  subgrain  microstructure  formation)  at  the 
mesoscale  and  (iii)  chemical  kinetics  at  the  molecular  scale.  The  centerpiece  of  the  model  is  an 
explicit  construction  that:  i)  gives  the  effective  or  macroscopic  behavior  of  plastically  deforming 
crystals  with  microstructure,  and  ii)  enables  the  reconstruction  of  optimal  microstructures  from 
the  computed  macroscopic  averages.  An  intrinsic  feature  of  the  optimal  deformation 
microstructures  is  the  presence  of  highly  localized  regions  of  plastic  deformation,  or  slip  lines 
(see  Fig.  2.1).  Temperatures,  strain-rates  and  pressures  in  these  slip  lines  rise  well  in  excess  of 
the  average  or  macroscopic  values.  Slip  lines  thus  provide  a  plentiful  supply  of  likely  initiation 
sites,  or  hot-spots,  in  defect- free  crystals. 

Relevance  to  Army 

This  project  is  part  of  an  integrated  research  program  to  address  various  fundamental  issues 
critical  to  the  development  of  insensitive  munitions  for  Future  Force  Insensitive  Munitions 
(FFIM).  The  developed  methodologies  should  provide  critical  data  on  atomistic  mechanisms  of 
initiation  and  detonation  in  condensed-phase  HE  and  support  the  development  of  advanced 
energetic  compositions.  Use  of  predictive  simulation  technique  should  decrease  the  time  and  cost 
of  development  for  enhancing  performance  and  safety  of  new  HE  and  formulations  for  FFIM. 
The  proposed  methodology  may  be  also  helpful  in  prediction  of  the  properties  and  responses  of 
HE  in  improvised  explosive  devices  (IED),  and  currently  is  being  used  by  Israeli  collaborator  at 
the  Hebrew  University  of  Jerusalem  (Prof.  Ronnie  Kosloff).  Besides  of  the  military  applications, 
developed  methodology  of  reactive  simulation  and  nanoscale  experimental  techniques  are 
already  being  used  to  study  reactive  processes  in  combustion,  catalysis,  fuel  cells,  and  other 
civilian  applications  of  the  reaction  modeling  in  nanostructures. 

Our  multiscale  simulations  of  hot  spots  formation  in  homogeneous  HE  due  to  the  stress 
localization  and  temperature  increase  at  the  slip  lines  in  subrain  microstructures  allows  to 


develop  a  predictive  model  of  initiation  and  sensitivity  of  realistic  HE  compositions  which  is 
impossible  to  construct  on  a  basis  of  conventional  mean-field  approaches  currently  used  for 
detonation  simulations  at  the  macroscale  level.  The  developed  codes  will  help  DoD  researchers 
in  the  the  development  of  new  and  improved  materials  and  structures  with  enhanced  energy 
density  and  reduced  sensitivity  for  a  wide  range  of  DoD  applications. 

Accomplishments 

1.  Large-scale  ReaxFF  simulations  of  the  shock  wave  propagation  in  a  single  crystal  of  PETN  to 
model  the  effect  of  initiation  of  chemical  reactions  and  shock-to-detonation  transition  in 
homogeneous  HE  (Figs.  1,  2,  4).  By  looking  at  particular  chemical  species  concentrations  as  a 
function  of  time  and  space  (Figs.  5-7),  we  have  been  able  to  correlate  the  observed  changes  in 
shock-front  velocity  with  the  storage  and  release  of  chemical  energy. 

2.  Predictive  modeling  of  anisotropy  of  shock  sensitivity  of  HE  single  crystals.  Compressive 
Shear  Reactive  Dynamics  (CSRD)  protocol  was  developed  to  study  the  effect  of  crystal  structure 
and  slip  directions  on  chemical  initiation  and  sensitivity  of  HE  single  crystal  (Figure  8).  We 
show  that  shocks  with  simultaneous  shear  and  compression  leads  to  enhanced  chemical  initiation 
and  that  in  some  materials  (PETN)  this  depends  dramatically  on  the  slip  system.  Our  protocol 
was  validated  by  comparison  with  experimental  results  of  J.  J.  Dick  on  single  crystal  of  PETN 
(Fig.  10  and  Table  1).  Now  the  protocol  is  being  used  to  simulate  other  single-crystal  HE  as 
RDX  (figures  1 1-15)  as  well  as  more  complex  HE  composites  with  grains  and  polymer  binders. 

3.  Empirical  vdW-correction  of  the  density  functional  theory  (DFT)  to  include  accurate  London 
Dispersion  for  calculations  of  wide-range  EOS  for  HE  crystals  with  accurate  description  of  both 
equilibrium  and  high-pressure  properties.  It  has  been  successfully  validated  on  hydrocarbon 
crystals  commonly  used  for  testing  purposes  (Figs.  1.1-1. 3).  We  now  use  it  to  determine  EOS  of 
HE  crystals  (PETN,  RDX,  HMX,  TNT)  and  provide  accurate  training  data  for  a  development  of 
new  generation  of  the  reactive  force  fields. 

4.  Fast  multiscale  models  of  shock-induced  plasticity  and  hot  spot  formation  in  energetic 
polycrystals.  We  have  assessed  the  feasibility  of  the  slip-line  initiation  mechanism  by  simulating 
a  PETN  plate  impact  experiment  and  comparing  the  resulting  predictions  with  experimental  pop- 
plot  data.  The  computed  characteristic  exponents  are  in  the  ballpark  of  experimental  observation, 
which  furnishes  a  modicum  of  validation  of  the  model  and  illustrates  the  ability  of  the  multiscale 
model  to  make  contact  with  full-scale  experimental  data  and  applications. 

We  have  specifically  assessed  the  ability  of  the  model  to  predict  experimentally  observed  pop- 
plot  characteristic  exponents  for  PETN.  We  exercise  the  multiscale  model  described  in  the 
foregoing  in  a  simple  plate  impact  configuration,  Fig.  2.3,  where  a  target  plate  of  polycrystalline 
PETN  that  is  impacted  on  its  rear  surface  by  a  rigid  striker  plate.  The  finite  element  mesh  used  in 
the  simulations  to  resolve  the  polycrystalline  structure  of  the  material,  is  shown  in  Fig.  2.3b.  The 
dimensions  of  the  sample  are  800  x  400  x  400  pm,  and  the  maximum  grain  size  is  100  pm.  Fig. 
2.5  shows  the  evolution  of  the  macroscopic  temperature  for  the  impact  velocity  of  700  m/s.  At 
material  points  where  deformation  microstructures  form,  the  macroscopic  temperature  follows  as 
the  volume-average  of  the  microscopic  temperature  field,  which  varies  sharply  on  the  scale  of 
the  microstructure.  A  typical  temporal  evolution  of  the  temperature  and  reacted  molar  fraction 
within  a  slip  line  is  shown  in  Fig.  2.7. 

The  comparison  with  experimental  data  is  predicated  on  the  assumption  that  initiation  in  a  pop- 
plot  test  requires  the  formation  of  a  critical  number  Nc  of  hot-spots.  Fig.  9  shows  the  pop-plot 
predicted  in  the  manner  just  described,  with  hot  spots  identified  with  slip  lines  that  attain  a 


temperature  greater  than  700  K,  a  reacted  fraction  greater  than  1%  and  a  pressure  greater  than  6 
GPa.  As  may  be  seen  from  the  Fig.  2.9,  the  pop-plot  predicted  by  the  model  is  consistent  with 
power-law  scaling,  in  agreement  with  experiments.  The  characteristic  exponent  predicted  by  the 
model  is  2.9,  which  slightly  overestimates  the  experimental  range  of  2.01  -  2.58  for  single 
crystal  PETN. 

Collaborations  and  Technology  Transfer 

1 .  The  Caltech  team  has  established  close  collaboration  with  personnel  from  the  Army  Research 
Laboratory  (ARL)  and  the  Naval  Research  Laboratory  (NRL).  Since  the  beginning  of  the  MURI, 
there  has  been  an  expanded  effort  to  exchange  research  data  and  codes  between  Caltech  team 
and  research  groups  at  ARL  and  NRL.  Caltech’s  simulation  codes  and  ReaxFF  force  fields  for 
reactive  atomistic  modeling  have  been  transferred  to  ARL,  NRL,  Air  Force  Institute  of 
Technology  at  Wright  Patterson  Air  Force  Base  (WPAFB),  Los  Alamos  National  Laboratory 
(LANL),  Sandia  National  Laboratories  (SNL),  NASA  Ames  Research  Center,  and  other  national 
research  institutions.  A  close  working  relationship  has  been  developed  with  the  High 
Performance  Technologies  Inc.  (HPTi)  at  WPAFB  on  customization,  distribution,  and  training 
for  ReaxFF  codes  and  force  fields  for  the  users  at  ARL,  NRL,  and  other  DoD  high-performance 
centers  (HPC).  The  ReaxFF  codes  and  force  fields  for  HE  are  being  transferred  to  DoD  labs  in 
collaboration  with  ARL  and  HPTi  researchers  (Betsy  Rice,  Anthony  Yau,  James  Lill,  James 
Larentzos).  Currently,  there  is  an  ongoing  HPTi/WPAFB  -  Caltech  effort  (PETTT  project#  PP- 
CCM-KY02-008)  on  the  support  and  further  modernization  of  ReaxFF  with  novel  optimization 
algorithms  developed  at  the  ARL  and  HPTi. 

2.  There  have  been  a  series  of  visits  to  ARL,  NRL,  HPTi,  LANL,  SNL,  other  universities,  and 
continuing  technical  interactions  on  the  multiscale  reactive  modeling  in  energetic  materials  with: 

-  Betsy  Rice  (ARL),  Jan  Andzelm  (ARL),  Anthony  Yau  (HPTi),  James  Lill  (HPTi),  James 
Larentzos  (HPTi):  optimization  and  training  of  the  ReaxFF  force  fields  for  energetic 
materials  and  applications; 

-  Peter  Chung  (ARL),  Jaroslaw  Knap  (ARL),  Betsy  Rice  (ARL):  coupling  continuum  (Prof. 
Ortiz)  and  atomistic  (Prof.  Goddard)  simulations  of  polycrystalline  plasticity,  hot  spots 
formation  and  initiation  of  detonation  in  heterogeneous  poly  crystalline  energetic  materials; 

-  Carter  White  (NRL),  Ivan  Oleynik  (USF):  anisotropic  nonlinear  elasticity  and  equations  of 
states  of  crystalline  EM  (PETN,  RDX,  HMX,  TATB,  nitromethane); 

-  James  Lill  (HPTi),  Betsy  Rice  (ARL),  and  Dr.  Aidan  Thompson  (SNL):  large-scale 
simulations  of  shock-induced  chemistry  in  HE  single  crystals  (RDX,  HMX.  PETN), 
polycrystalline  HE  and  HE/binder  composites  (RDX/HTPB,  PETN/HTPB)  using  DoD 
HPCMP  high-performance  computers  (CAP  project); 

-  Ed  Kober  (LANL)  and  Tomas  Sewell  (Univ.  of  Missouri):  ReaxFF  reactive  modeling  of 
chemical  kinetics  of  thermal  decomposition  and  formation  of  hot  spots  in  HMX,  TATB,  and 
nitromethane; 

-  Aidan  Thompson  (SNL):  shock-to-detonation  transition  in  PETN  and  CL-20; 

-  Ronnie  Kosloff  (Hebrew  University  of  Jerusalem):  ReaxFF  modeling  of  the  HE  response  to 
the  terahertz  (THz)  spectroscopy  probing;  simulation  of  initiation  and  decomposition  kinetics 
in  nitromethane,  TNT,  TATP,  and  DATP. 


3.  Prof.  William  Goddard  gave  multiple  invited  talks  at  national  and  international  workshops  and 

symposiums  on  reactive  force  fields  and  RMD  modeling  of  energetic  materials  and  sensitivity: 

-  Gordon  Research  Conference  on  Energetic  Materials,  Tilton,  NH,  June,  2006. 

-  7th  International  Conferences  “New  Models  and  Hydrocodes  for  Shock  Wave  Processes  in 
Condensed  Matter”  2008  (Lisbon,  Portugal) 

-  Goddard  all  day  visit  and  seminar  LANL,  April  19,  2010 

-  Goddard  all  day  visit  and  seminar  SNL,  April  20,  2010 

-  International  Conference  and  Advanced  School  Turbulent  Mixing  and  Beyond  TMB-2009, 
Trieste  Italy,  June  2009,  Goddard  gave  invited  talk 

-  Prof.  Michael  Ortiz  visited  ARL  multiple  times  (2009,  2010)  to  give  a  presentation  on 
multiscale  modeling  of  HE  detonation  initiation  with  the  finite  elements  code  developed  at 
Caltech.  The  code  is  being  developed  with  a  goal  to  transfer  it  to  ARL  for  simulations  of 
impact  sensitivity  and  initiation  in  polycrystalline  and  polymer-bonded  explosives. 

-  Dr.  Sergey  Zybin  and  Dr.  Adri  van  Duin  visited  ARL  (2005,  2007)  and  SNL  (2006,  2008)  to 
discuss  development,  maintenance,  and  use  of  the  ReaxFF  force  fields  in  EM  applications  as 
well  as  parallel  implementation  of  ReaxFF  in  multiprocessor  software  Grasp  and  Lammps. 
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Figure  1.  Tetragonal  PETN  unit  cell  with  a  schematic  of 
shock  generation  by  impact  against  a  stationary  wall  at 
particle  (i.e.  piston)  velocity  -up 


time  for  a  weak  shock  (3  km/s)  and  a  strong  shock  (4  km/s)  along  [100] 
direction  in  PETN  crystal  with  density  ~1.5  g/cm3.  The  observed 
acceleration  of  the  shock-wave  (shock  to  detonation  transition)  is  due  to 
the  net  energy  release  that  occurs  when  final  reaction  products  are 
formed.  In  the  case  of  the  weak  shock,  the  initial  deceleration  is  due  to 
energy  storage  (endothermicity)  in  intermediate  reaction  products. 


Figure  2.  Hugoniot  adiabat  from  ReaxFF-MD  simulations  of 
unidirectional  shocks  in  PETN  single  crystals  with  density 
-1.7  g/cm3  show  good  agreement  with  experiment  [Error! 
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Figure  5.  Distribution  of  primary  reaction  products  (dissoci¬ 
ated  N02  per  PETN  molecule)  along  [100]  shock  direction. 
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Figure  6.  Distribution  of  secondary  reaction  products  (H20 
per  PETN  molecule)  along  the  [100]  shock  direction. 


Figure  7.  Distribution  of  secondary  reaction  products  (N2  per 
PETN  molecule)  along  the  [100]  shock  direction. 


S  ch.eiTia.tic  s  o  f  D  simulation  of  constant -rate  hear 
along  one  of  the  slip  plane  in  FEIN  single  crystal  to  evaluate 
the  effect  of  hear  deformation  on  chemical  initiation. 


plane 


Shear  simulations  on  10%-compressed  PETN  (after  6ps) 


Table  1 
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I 
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0.2 

0.49 

0.69 
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a>  Shock  plane  normal  is  parallel  to  the  shock  direction. 

b>  S  -  sensitive,  I  -  insensitive,  M  -  intermediate  and  velocity-dependent. 

c>  T0  is  the  initial,  Tm  is  the  maximum,  and  $c  is  the  steady  state  shear  stress 


Insensitive 


ye  nsitive 


:v-ID  simulation  of  constant  -rate  hear 
predicts  much  faster  chemical  initiation  (dissociation  of  NO;) 
in  PETN  single  crystal  for  experimentally  more  sensitive 
[1 10]  and  [D01  ]  directions  of  hock  compression. 


Least  Hindered  Shear  System  (B)  in  RDX 


Least  Hindered  Shear  System  (B)  in  RDX 


Fig-11 


Steric  hindrance  index:  0.09 


Fig.12 


Steric  hindrance  index:  0.09 
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Constant-rate  shear  on 
20%-compressed  RDX 
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normal  to  shear  plane. 
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I  Shear  deformation  sustains  compressed  pressure 
I  1  negative;  along  hindered  shear  direction,  while 
■  builds  up  expanded  pressure  (positive)  along  direction 
I  normal  to  shear  plane. 
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Fig.1.1  Dispersion  correction  (E±j^)  for  a  pair  of  carbon  atoms 
developed  in  t£ its  work  (denoted  as  LG)  compared  with  those  in 
literatures  including  WYl“  WY2M,  EHFSK".  and  OBS:&  The 
same  C-C  equilibrium  bond  distance  Ro  =  3.851  A  and  Q  =  5Sd.8 
kcalsnol  A*  are  used  for  die  sake  of  comparison.  C^'r 

without  damnum  is  also  nlotted  as  a  reference 
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Fig.  1.3  (a)  EDS  of  benzene  crystal  (orthorhombic  phase  I) 
and  (b)  naphthalene  crystal  calculated  using  PBE  and  PBE-D 
together  with  the  extrapolated  0  K  EOS.  The  inserts  show  the 
unit  cells  of  the  crystals. 
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I  Fig.1.2  Binding  energy  of  benzene  dimers  as  a  fimction  of  distance  calculated  using  pure  PBE  and  its  dispersion  corrections  (PBE- 
D)  for  (a)  T-shaped,  (b)  sandwich,  (c)  parallel-displaced  configurations. 


FIG.  1:  Simplest  possible  microstructure.  Macroscopic  defor¬ 
mation  i3  has  two  active  slip  systems .  and  decomposes  into 
first  order  laminates  3i  and  ih  with  single  and  double  slip, 
respectively.  Similarly.  decomposes  into  second  order  lam¬ 
inates  p'ife  and  3i  with  single  and  double  shp,  respectively 
Highly  localised  plastic  deformations  occur  in  tii  which  isccai- 
sidered  as  a  potential  hot  spot. 
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FIG.  3 1  Plate  impact  test  simulation,  a)  Definition  of  computational  domain  and  boundary  conditions,  b)  Finite-element  mesh 
resolving  the  poly  crystalline  structure  of  the  sample. 
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Fig. 2. 5  Macroscopic  evolution  of  temperature  for  an  impact  velocity  of  rQOm/s.  Red  and  blue  colors  in  the  figure  correspond 
to  temperatures  of  £?=40D  K  and  S=273  K,  respectively. 
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.  ,  ...  Fig.2.9:  PEI  N  plate  impact  simulation.  Power-law  scaling 

Fig.2.7:  temperature  and  reacted  fraction  evolution  m  a  typ-  i  .■  u  ,  .i  «-  i  r  .  ,  i  (1 

,  .  .  ^  ,  L  r  relation  between  the  inverse  or  number  oi  hot  spots  and  the 

ical  hot-  spot.  Red  diamonds  and  blue  circles  correspond  to-  .  i  ■.  Ti  ■  i  .i  ■  i  ^  t 

r  /  impact  velocity.  1  he  circles  are  the  simulation  results  lor 

the  temperature  and  the  reacted  fraction^  respect  ivelv.  i-jt  ■  .  i  -  ■  ti  i  *.  c**4.uj* 

r  r  dine  rent  impact  velocities,  ihe  least  square  ht  to  the  data 

shown  by  the  solid  Hne  has  the  slope  of  2,91. 


